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OPTIMAL  STAGING  AND  SCHEDULING 
IN  AIRLIFT  OPERATIONS 


INTRODUCTION 

The  staging  problem  in  airlift  operations  is  a  generalized  resource  allo¬ 
cation  problem  of  achieving  maximum  aircraft  utilization  by  optimally  staging 
aircrews  at  different  airbases  at  the  start  of  an  airlift.  The  problem  is 
inherently  dynamic  in  nature  and  involves  time  delays.  In  this  paper  we  will 
demonstrate  how  to  model  a  simplified  version  of  this  problem  as  a  time- 
expanded  graph  [1]  and  solve  the  resulting  mixed-integer  programming  problem. 
The  basic  approach  employed  here  is  a  generalization  of  that  presented  by  Ford 
and  Fulkerson  [1]  and  Gearhart  [2].  In  particular,  we  address  the  problem  of 
multiple  routes,  each  of  which  may  form  a  nonsimple  loop.  We  also  address  the 
related  problem  of  ordering,  or  prioritizing,  routes  in  a  linear  programming 
envi ronment. 


THE  PROBLEM 

The  basic  system  under  consideration  consists  of  P  aircraft,  C  aircrews, 
and  B  airbases:  B ( 1 ) ,  B ( 2 ) ,  ....  B(B).  B(l)  is  the  home  base.  This  is  where 
all  the  P  aircraft  are  initially  located  and  where  all  routes  start  and  end. 
A  route  is  any  finite  sequence  of  airbases 

hj ) ,  B(n2),  . . . ,  B(nk). 

Here  any  consecutive  pair  (B(n^),B(n.+  ))  is  called  a  leg  of  the  route.  The 

only  restrictions  on  any  route  are  that  the  first  and  last  stops  should  be  the 
home  base  and  that  no  other  stop  should  be  there.  Note  that  not  all  airbases 
need  to  be  included  in  a  route  and  that  an  airbase  may  appear  more  than  once 
in  the  same  route. 

A  mission  is  a  route  with  the  associated  attributes  of  aircraft  and 
starting  time.  The  number  of  routes  in  our  system  is  fixed,  but  the  number  of 
missions  can  vary  since  more  than  one  mission  can  follow  the  same  route.  The 
number  of  missions  depends  on  the  duration  of  the  airlift  operations  under 
consideration. 

Each  mission  is  flown  by  a  single  aircraft;  however,  it  is  generally 
flown  by  more  than  one  aircrew,  each  responsible  for  one  leg  of  the  mission. 
An  aircrew  is  required  to  rest  for  a  certain  period  of  time  after  each  leg  of 
a  mission,  but  the  mission  may  be  continued  whenever  any  rested  aircrew  is 
available.  If  a  rested  aircrew  is  unavailable  at  an  airbase  B(n)  to  replace 
another  (incoming)  aircrew  at  the  end  of  the  leg  (B(m),B(n)),  the  aircraft, 
and  hence  the  mission  in  question,  will  have  to  be  delayed  until  a  rested  air¬ 
crew  does  become  available.  In  this  simplified  model,  we  have  chosen  to 
ignore  a  few  restrictions  both  on  the  aircraft  and  the  aircrews;  some  are 
ignored  to  keep  the  model  simple,  and  others  because  their  inclusion  would 
render  the  problem  unsolvable  by  reasonable  analytical  methods.  Typical  of 
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the  latter  is  the  limitation  on  the  total  flying  hours  of  the  aircrews  during 
any  30-day  period;  a  model  including  this  restriction  would  have  to  include  a 
running  total  of  the  flying  hours  of  each  aircrew.  For  a  simulation  model, 
this  can  be  handled  quite  easily  (for  example,  using  SIMSCRIPT,  although  not 
Q-GERT  nor  GPSS);  but  for  an  analytic  optimization  model,  this  restriction  is 
impractical  to  include. 

The  problem  that  we  are  considering  here  can  be  stated  as  follows:  Given 
C  aircrews,  P  aircraft,  R  routes,  an  ordered  sequence  of  missions  (to  be  elab¬ 
orated  later),  and  the  duration  of  the  airlift  (time  period  T),  how  should  the 
aircrews  be  initially  distributed  among  the  bases  so  that  the  flying  hours  by 
all  aircraft  are  maximized  during  the  airlift?  The  ordered  sequence  of 
missions  alluded  to  is  implicitly  given  by  specifying  the  proportion,  p(i) 
(i=l,  ...,R),  of  the  use  of  each  route. 


THE  TIME-EXPANDED  GRAPH  FORMULATION 

The  basic  idea  as  suggested  by  Gearhart  [2]  is  to  model  the  system  as  a 
time-expanded  graph  [1].  By  assuming  that  all  changes  in  the  system  take 
place  at  integer  multiples  of  a  basic  time  unit,  which  we  will  take  as  1  for 
simplicity,  we  can  discretize  the  time  period  (0,T)  so  that  the  system  is 
fully  described  by  its  states  at  times  0,  1,  2,  ....  T.  Here  T  can  be  assumed 
to  be  an  integer.  As  we  shall  see,  our  system  can  be  represented  by  the  graph 
G=(N,A),  where  N  is  a  finite  set  of  nodes  and  A  is  a  finite  set  of  arcs  join¬ 
ing  pairs  of  nodes.  The  graph  is  completely  specified  if  we  define  all  the 
nodes,  all  the  arcs,  all  the  flows  along  the  arcs,  and  all  the  constraints. 
We  will  proceed  to  do  this. 


The  Nodes 

Four  types  of  nodes  comprise  N:  route,  aircrew,  schedule,  and  source 
nodes.  We  first  define  the  set  of  route  nodes  associated  with  time  n  as  R(n). 
The  route  nodes  represent  airbases  B(l),  B(2),  ...,  B(B),  and  possibly  copies 
of  them— depending  on  the  route  structure  in  the  system.  Specifically, 
suppose  the  system  has  R  routes  and  that  the  i-th  route  has  L(i)  legs;  then 
R(n)  contains  exactly  L(l)  +  L(2)  +  ...  +  L(R)  -  R  +  1  route  nodes,  each  of 
which  (with  one  exception)  corresponds  to  a  unique  stop  in  a  unique  route. 
The  same  airbase  may  have  to  be  represented  by  more  than  one  node  to  ensure 
that  aircraft  will  adhere  to  their  routes.  If  airbase  B(i)  is  the  j-th  stop 
in  the  r-th  route,  the  node  representing  B(i)  at  time  n  will  be  B(i,n,r,j). 
The  exception  to  a  route  node  being  a  unique  stop  in  a  unique  route  is  the 
node  for  home  base,  B(l).  At  time  n,  the  common  node  for  B(l)  will  be  B(l,n), 
and  home  base  will  be  identified  by  B(l,n,k,l)  and  by  B(l,n,k,L(k)+l)  ( k=l , 
...,R).  The  set  of  route  nodes  is  the  union  of  R(l),  R(2),  ...»  and  R(T). 

The  family  of  aircrew  nodes  consists  of  B  nodes  at  each  time  n,  each 
representing  the  aircrew  pool  at  a  specific  airbase.  These  nodes  will  be 
denoted  by  C { i , n )  (i  =1,...,B;  n=0,l, ...,T). 

The  family  of  schedule  nodes  consists  of  a  master  schedule  node,  M,  and  R 
subsidiary  schedule  nodes,  M(i,n)  (i=l,...,R),  for  each  time  n  (n=0,...,T). 
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The  master  schedule  node  is  used  to  control  the  order  in  which  the  different 
mission  types  are  generated,  and  the  subsidiary  schedule  nodes  are  used  as 
counters  for  the  different  types  of  missions  that  have  been  started  on  or 
before  time  n. 

Finally  the  family  of  source  nodes  consists  of  exactly  one  node,  S, 
representing  the  number  of  aircraft  and  aircrews  in  the  system. 


The  Arcs 

Arcs  are  ordered  pairs  of  nodes.  An  arc  is  of  interest  only  if  a  flow 
(along  the  arc)  representing  a  meaningful  quantity  exists.  Four  types  of  arcs 
are  in  our  system:  arcs  between  route  nodes,  between  aircrew  nodes,  between 
schedule  nodes,  and  between  source  and  other  nodes.  We  will  elaborate  on 
these  now. 

Arcs  Between  Route  Nodes--Three  types  of  arcs  connect  the  route  nodes. 
These  correspond  to  the  three  states  an  aircraft  can  be  in--maintenance, 
delayed,  and  flying.  The  first  type  is  of  the  form 

(B(l ,n) ,B(l,n+m)) , 

where  in  is  the  aircraft  maintenance  time  at  home  base.  The  flow  M(n)  on  one 
of  these  arcs  represents  the  number  of  aircraft  starting  maintenance  at  time 
n.  Recall  that  B(l,n,k,l)  and  B(l,n,k,L(k)+l)  (k=l,...,R)  both  represent  the 
same  node  (home  base  at  time  n)  and  are  represented  by  8(1, n). 

The  second  type  of  arc  connecting  the  route  nodes  is  of  the  form 

(B(i ,j,k,h) .B(i ,j+l,k,h)). 

The  flow  D(i,n,k,h)  on  one  of  these  arcs  represents  the  number  of  aircraft 
being  delayed  (for  example,  waiting  for  a  rested  aircrew)  at  airbase  B(i)  at 
time  n;  the  aircraft  in  question  being  at  the  h-th  stop  on  the  k-th  route.  In 
the  special  case  in  which  i  =  1  (home  base),  D(l,n,k,l)  (k=l,...,R)  represents 
e  same  flow  and  is  denoted  by  D(l,n). 

Ten  T  -'st  type  of  arc  connecting  the  route  nodes  is  of  the  form 

(B(i,n,k,h),B(b(k,h+l),n+w(k,h),k,h+l)); 

here  b(k,h)  is  the  h-th  stop  on  the  k-th  route,  and  w(k,h)  is  the  flight  time 
(an  integer)  between  b(k,h)  and  b(k,h+l).  The  flow  F(i,n,k,h)  on  one  of  these 
arcs  represents  the  number  of  aircraft  flying  from  b(k,h)  to  b(k,h+l),  start¬ 
ing  at  time  n  and  arriving  at  b(k,h+l)  at  time  n  +  w(k,h). 

Arcs  Between  Aircrew  Nodes--Three  types  of  arcs  connect  the  aircrew 
nodes.  These  correspond  to  the  three  states  an  aircrew  can  be  in--rest, 
delayed,  and  flying.  The  first  type  is  of  the  form 

(C(i,n),C(i,n*l)). 
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The  flow  d(i,n)  on  one  of  these  arcs  represents  the  number  of  aircrew  being 
delayed  (after  being  rested)  for  one  unit  of  time  at  airbase  B(i)  starting  at 
time  n. 

The  second  type  of  arc  connecting  the  aircrew  nodes  is  of  the  form 
(C(i ,n) ,C(i ,n+r) ) . 

The  flow  r(i,n)  on  ono  of  these  arcs  represents  the  number  of  aircrew  starting 
to  rest  (for  r  time  units)  at  time  n  on  airbase  B(i). 

The  last  type  of  arc  connecting  the  aircrew  nodes  is  of  the  form 

(C(i ,n) ,C(k ,n+t(i ,k) )) . 

Here  i  is  not  equal  to  k,  and  t(i,k)  represents  the  flight  time  between  B(i) 
and  B(k).  We  are  assuming  that  all  aircraft  are  of  the  same  type,  so  the 
flight  time  between  any  two  airbases  can  be  assumed  to  be  the  same,  indepen¬ 
dent  of  aircraft.  The  flow  f(n,i,k)  on  one  of  these  arcs  represents  the  num¬ 
ber  of  aircrews  flying  from  airbase  B(i)  to  B(k),  starting  at  time  n  and 
arriving  at  B(k)  at  time  n  +  t(i,k). 

Arcs  Between  Schedule  Nodes — Two  types  of  arcs  connect  the  schedule 
nodes”  The  first  type  is  of  the  form 

(M,M(i  ,n)) . 

These  are  arcs  between  the  master  schedule  node  and  the  subsidiary  schedule 
nodes.  The  flow  s(i,n)  on  one  of  these  arcs  represents  the  number  of  missions 
to  be  sent  out  on  route  i  at  time  n.  In  our  model,  we  assume  s(i,n)  is  either 
one  or  zero. 

The  second  type  of  arc  connecting  the  schedule  nodes  is  of  the  form 
(M(i ,n) ,M(i ,n+l) ) . 

As  mentioned  before,  the  subsidiary  schedule  nodes  serve  as  counters  for  the 
flown  missions  of  a  given  type.  The  flow  S(i,n)  on  one  of  these  arcs  repre¬ 
sents  the  total  number  of  missions  on  route  i  that  have  been  started  on  or 
before  time  N. 

Arcs  Between  the  Source  Node  and  Other  Nodes — Two  types  of  arcs  involve 
the  source  node,  S.  The  first  type  is  of  the  form 

(S,B(1,0) ) . 

The  flow  P  along  this  arc  represents  the  number  of  aircraft  in  the  system.  We 
assume  all  aircraft  had  their  initial  maintenance. 

The  second  type  of  arc  involving  the  source  node  is  of  the  form 

(S,C(i ,0)) . 
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The  flow  c(i)  on  one  of  these  arcs  represents  the  number  of  aircrews  initially 
staged  at  the  airbase  B(i).  We  assume  all  aircrews  are  rested  initially. 


The  Constraints 

There  are  three  types  of  constraints:  integer,  aircraft-aircrew,  and 
scheduling.  The  integer  constraints  force  the  flows  (variables)  c(i),  F(i,n, 
k,h),  and  s(i,n)  to  be  integers.  The  ai rcraft-ai rcrew  constraints  force  the 
number  of  aircrews  flying  out  of  any  airbase  to  be  equal  to  the  number  of  air¬ 
craft  flying  out  of  that  airbase.  Finally,  the  schedule  constraints  force  the 
right  missions  to  be  generated;  a  detailed  discussion  of  this  is  given  in 
Appendix  A.  One  of  the  assumptions  we  are  making  on  scheduling  is  that  at  any 
time  n,  the  number  s(n)  of  missions  generated  (regardless  of  type)  is  at  most 
one.  This  assumption  is  made  for  purely  technical  reasons;  without  it  we 
would  have  to  deal  with  an  increased  number  of  equations.  It  does  not  affect 
practical  application  since  the  time  unit  can  be  made  small  enough  for  this 
assumption  to  be  realistic. 


THt  MIXED-INTEGER  PROGRAMMING  FORMULATION 

By  balancing  the  flows  at  the  nodes  in  our  graph  (except  for  the  source 
node,  S,  and  master  node,  M)  and  taking  into  consideration  the  constraints  in 
the  system,  we  can  set  up  a  mixed- integer  programming  model  for  our  problem. 
To  be  more  explicit,  we  will  demonstrate  this  by  working  out  an  example.  We 
assume  there  are  three  airbases  (B=3)--B(l),  B(2),  and  B (3 ) -- and  two  routes 

(R-2): 


Route  1:  B(l),  B(2),  B(3),  B(2),  B(l) 

Route  2:  B(l),  B(3),  B(2),  B(3),  B(l). 

We  assume  the  flight  time  between  airbase  B(i)  and  B(j)  is  t(i,j),  air¬ 
crew  rest  time  is  r,  and  aircraft  maintenance  time  is  m.  We  also  assume  the 
actual  scheduling  starts  at  time  n  =  0,  so  c(i)  =  d(i,0)  (i=l,2,3).  Then  the 
equations  and  constraints  describing  this  model  are  given  as  follows.  (A  list 
of  the  variables  used  and  their  definitions  can  be  found  in  Appendix  B.) 

Source  node: 


d(l ,0)  +  d(2,0)  +  d(3 ,0)  =  C 
0(1,0)  =  P. 

Route  nodes  (for  n=l,2,...T): 

M(n)  -  F (2 ,n-t (2, 1 ) , 1 ,4)  +  F(3,n-t(3,l) ,2,4) 

D(l,n)  +  F ( 1 ,n,  l ,  1 )  +  F(l,n,2,l)  =  D(n-l)  +  M(n-m) 
D(2,n,l,2)  +  F(2,n,l,2)  =  D(2,n-l,l,2)  +  F (1 ,n-t(l ,2)  ,1,1) 

D(3,n,l,3)  +  F(3,n,l,3)  =  D(3,n-l,l,3)  +  F (2,n-t(2,3)  ,1,2) 

D(2,n,l,4)  +  F(2,n,l,4)  =  D(2,n-l,l,4)  +  F(3,n-t(3,2) ,1,3) 

D(3,n,2,2)  +  F(3,n,2,2)  =  D(3,n-1,2,2)  +  F(l,n-t(l,3)  ,2,1) 

D(2,n,2,3)  +  F (2,n,2,3)  =  D(2,n-1,2,3)  +  F(3,n-t(3,2) ,2,2) 

D(3,n,2,4)  +  F(3,n,2,4)  =  D(3,n-1,2,4)  +  F(2,n-t(2,3)  ,2,3) . 
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Aircrew  nodes  (for  n=l ,2, . . . ,T) : 

f (n,l  ,3)  +  f(n,l ,2)  +  d(l,n)  =  d(l,n-l)  +  r(l,n-r) 
r(l,n)  =  f(n-t(2,l)  ,2,1)  +  f(n-t(3,l)  ,3,1) 
f (n,2,l)  +  f(n,2,3)  +  d(2,n)  =  d(2,n-l)  +  r(2,n-r) 
r(2,n)  =  f ( n-t ( 1 ,2)  ,1,2)  +  f (n-t(3,2)  ,3,2) 
f(n,3,l)  +  f(n,3,2)  +  d(3,n)  =  d(3,n-l)  +  r(3,n-r) 
r(3,n)  =  f (n-t(l  ,3)  ,1,3)  +  f (n-t (2 ,3)  ,2 ,3) . 

Schedule  nodes  (for  n=l ,2, . . .  ,T) : 

S(n,l)-s(n,l)  =  S(n-1 ,1 ) 

S(n,2)-s(n,2)  =  S(n-1,2). 

Ai rcraft-ai rcrew  constraints  (for  n=l,...,T): 

f(n,l,2)  =  F( 1  ,n, 1 , 1 ) 

f (n ,1 ,3)  =  F ( 1 ,n ,2 , 1 ) 

f(n,2,l)  =  F(2,n,2,3) 

f(n,2,3)  =  F (2 ,n ,1 ,2)  +  F(2,n,2,3) 

f (n,3,l)  =  F(3,n,2,4) 

f (n,3,2)  =  F(3 ,n ,1 ,3)  +  F(3,n,2,2). 

Schedule  constraints  (for  n=l,...,T,  and  W  a  sufficiently  large  number): 

F(1  ,n,l  ,1)  =  s(n,l) 

F(1 »n,2 ,1 )  =  s(n,2) 

2p(2)S(n,l)  -  2p(l)S(n,2)  +  Ws(n,l)  <  W  +  p(l)  -  p(2) 

2p(l)S(n,2)  -  2p(2)S(n,l)  +  Ws(n,2)  <  W  +  p(2)  -  p(l) 
s(n)  -  s(n,l)  -  s(n,2)  =  0. 

Integer  constraints:  The  integer  variables  are 

d( 1 ,0) ,  d(2,0),  and  d(3,0); 

and  the  0-1  integer  variables  are 

F(1  ,n,l  ,1) ,  F(2  ,n  ,1 ,2) ,  F(3,n,l,3),  F(2,n,l,4),  F(l,n,2,l), 

F(3,n,2,2) ,  F(2,n,2,3) ,  F(3,n,2,4),  and  s(n). 

The  objective  function  of  this  mixed-integer  problem  is  the  total  flying 
hours  by  all  aircraft  and  is  given  by 

U  =  l  |F(l,n,l,l)  +  F(2 ,n ,1 ,2)  +  F(3,n,l,3)  +  F(2,n,l,4) 
n=l 

+  F(l,n,2,l)  +  F(3,n,2,2)  +  F(l,n,2,3)  +  F(3,n,2,4)). 

Thus  the  problem  of  optimizing  aircraft  utilization  reduces  to  maximizing  U 
subject  to  all  the  constraints  stated  in  this  section. 


RESULTS 

The  mixed-integer  problem  was  solved  on  the  UNIVAC  1108,  using  the  mixed- 
integer  programming  (MIP)  program  in  the  Functional  Mathematical  Programming 
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System  (FMPS)  [3].  This  computer  program  was  developed  by  Bonner  and  Moore 
Associates,  Inc.  for  Sperry  Rand  Corporation  and  is  based  on  the  ideas  of 
Forrest  et  al .  and  Mitra  [4,5],  Table  1  summarizes  the  results  of  several 
runs  using  different  sets  of  parameters.  In  these  runs,  T  =  12,  P  =  3,  p(l) 
*  0.3,  p(2)  =  0.7,  t(l,2)  =  1.0,  t( 1 ,3)  =  2.0,  t(2,3)  =  1.0,  and  m  =  3.0.  The 
only  parameters  that  we  changed  from  run  to  run  were  C,  the  number  of  air¬ 
crews,  and  r,  the  aircrew  rest  time.  The  solution  to  the  mixed-integer  prob¬ 
lems  includes  aircraft  utilization,  U  (defined  here  as  the  total  number  of 
basic  time  units  flown  by  all  aircraft  during  T  basic  time  units),  and  the 
initial  distribution  of  aircrews:  d(l,0),  d(2,0),  and  d(3,0)  at  airbases 
B(l),  B(2),  and  B(3)  respectively.  Recall  that  the  number  of  aircrews  ini¬ 
tially  staged  at  an  airbase  is  represented  by  c(i)  =  d(i,0)  (i  =  l,2,3).  (The 
quantity  U  shown  in  Table  1  is  a  slight  overestimate  of  the  actual  aircraft 
utilization,  since  we  have  credited  the  complete  flight  time  of  a  leg  to  an 
aircraft  whenever  it  starts  that  leg;  this  induces  overestimates  on  the  flight 
times  of  all  aircraft  that  are  still  flying  at  time  T.)  The  dimension  of  the 
matrix  that  appeared  in  the  MIP  program  in  all  cases  is  330  rows  (constraints) 
and  409  columns  (variables),  117  of  the  latter  being  integers.  The  number  of 
iterations  (Iter)  and  CPU  time  shown  in  the  table  are  dependent  on  the  so- 
called  cutoff  value  for  that  run;  it  is  essentially  a  lower  bound  (provided 
by  the  user)  for  the  objective  function.  A  wisely  chosen  cutoff  value  can 
reduce  dramatically  the  number  of  nodes  to  be  searched  in  the  branch-and-bound 
procedure,  and  hence  the  number  of  iterations  and  CPU  time. 


TABLE  1.  RESULTS  OF  OPTIMIZATION 


RUN 

c 

r 

U  d 

(1,0) 

d(2 ,0 )  d 

1(3,0) 

Iter 

CPU 

(min) 

1 

4 

2 

16 

2 

0 

2 

1195 

9.98 

2 

4 

3 

★ 

★ 

* 

★ 

* . 

★ 

3 

4 

4 

14 

3 

0 

1 

2758 

16.79 

4 

5 

2 

24 

3 

1 

1 

1136 

9.16 

5 

5 

3 

20 

3 

1 

1 

1757 

15.41 

6 

5 

4 

14 

4 

0 

1 

3150 

18.98 

7 

6 

2 

24 

3 

1 

2 

2796 

23.66 

3 

6 

3 

21 

3 

1 

2 

4213 

38.46 

9 

6 

4 

18 

3 

1 

2 

3008 

20.98 

r — 
U— 

No.  of  aircrews  in  the  system. 

Aircrew  rest  time. 

Total  hours  flown  by  all  aircraft  during  the 

ai  rl  i  ft . 

d(l ,0) , 

d(2,0). 

d  ( 3 , 0 )  - 

-Initial  distribution  of 

ai rcrews 

at  airbases. 

Iter--No.  of  iterations. 

*The  Functional  Mathematical 

Programni  ng 

System 

[3]  fails  to  yield 

sol ut i on. 
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CONCLUSION 

In  this  paper  we  have  demonstrated  the  feasibility  of  using  mixed- integer 
programming  techniques  to  estimate  the  optimal  staging  policy  in  Air  Force 
airlift  operations.  At  present  the  method  suffers  from  the  dimensionality  of 
the  problem;  the  number  of  variables  and  constraints  increases  rapidly  with 
the  number  of  aircrews,  the  number  of  airbases,  the  number  of  routes,  and  most 
importantly,  the  duration  of  airlift  operations.  In  an  actual  airlift  prob¬ 
lem,  each  of  these  parameters  is  generally  tenfold  larger  than  the  one  we 
have  considered  here.  However,  due  to  the  sparsity  and  the  structure  of  the 
matrix  involved,  there  may  be  powerful  (decomposition)  techniques  that  can  be 
used  to  circumvent  this  difficulty.  This  is  an  area  for  further  research. 
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AJ’I’LNDI X  A:  MISSION  SCHEDULING  (ORDERING) 


In  this  appendix,  we  will  show  how  the  proportion,  p(i),  of  route  usage 
can  lead  naturally  to  mission  scheduling  (ordering).  Here  we  are  given  a 
finite  sequence  of  proportion  p(i)  (i=l,...,R)  so  that 

0  <  p(i )  _<  =  1 ,  and 

p(l)  +  p(2)  +  ...  +  p(R)  =  1 

where  R  is  the  number  of  distinct  routes*  Let  S(n)  be  the  total  number  of 
missions  that  have  been  sent  out  before  time  n+1,  and  similarly  let  S(n,k)  be 
the  total  number  of  missions  of  type  k  (k=l,...,R)  that  have  been  sent  out 
before  time  n+1.  Thus 

S(n)  -  S (n ,  1 )  +  S(n,2)  +  ...  +  S(n,R).  (Al) 

Ideally,  we  would  like  to  have  S(n,k)  =  p(k)S(n)  (k=l,...,R)  for  all  n.  But 
since  S(n,k)  and  S(n)  are  integers,  these  conditions  are  seldom  satisfied.  We 
can  view  S(n,k)/p(k)  (k~l,...,R)  as  R  different  estimates  of  S(n).  Our  prob¬ 
lem  is  then  to  find  a  method  for  selecting  the  sequence  of  missions  to  be  sent 
out  so  as  to  minimize  the  deviations  ($(n,k)/p(k)-S(n) )  (k=l,...,R)  at  each 
time  n. 

Note  that  from  equation  (Al)  we  have 

S(n)  =  (p(l)+  — +p(R) ) 

$(n)  -  S(n,k)  +  _  +  S(n,R) 

cr 

p( 1 ) (S(n, l)/p( 1 )-S(n) )  +  ...  +  p(R)(S(n,R)/p(R)-S(n))  =  0. 

That,  is,  the  average  deviation  is  zero  for  all  n.  A  criterion  to  use  to 
insure  small  deviations  (from  S ( n ) )  is  that  the  average  absolute  deviations 
from  S(n)  be  as  small  as  possible  for  each  n.  Another  criterion  is  that  the 
average  squared  deviation  should  be  minimized. 


Minimization  of  Average  Absolute  Deviation 

Assuming  S(n,k)  (k=l,..*,R)  has  already  been  specified,  the  selection  of 
the  next  type  of  mission  should  be  based  on  minimizing  the  average  absolute 
deviation  (from  S(n^l))  at  time  n+1;  i.e., 

I  (l)}S(n+l,l)/p(l)  -  S (n+1) I  +  p(2)  |s(n+l,2)./p(2)  -  S (n+1)  j 

+  ...  +  p(R)|S(n+ l,R)/p(R)  -  S(n+1) | .  (A2) 

Note  that  if  route  i  is  picked  at  time  n+1,  then 

S ( n f 1 , k )  S(n,k)  +  v ( i , k  ) 
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and 


S(n+1)  =  S(n)  +  1. 

Here  v(i,k)  =  1  if  i  =  k,  and  0  otherwise.  Incidentally,  S(n)  need  not  be  n 
in  general,  since  if  an  aircrew  or  an  aircraft  is  not  available,  no  mission 
can  be  sent  out,  and  in  that  case  S(n+1)  =  S(n)  and  S(n+l,k)  =  S(n,k)  (k-1. 
The  following  lemma  gives  the  selection  rule  based  on  minimizing 
expression  (A2). 

Lemma  A .  Let  Z(n+l,k)  (k=l,...,R)  be  the  average  absolute  deviation  from 
S(n+1 )--i ,e. ,  the  value  of  expression  (A2)  —  if  a  mission  using  the  k-th  route 
is  chosen  at  time  n  +  1.  And  let  e(n+l,k)  =  $(n,k)  -  p(k)S(n+l);  this  is 
related  to  the  difference  between  the  actual  and  the  expected  number  of  mis¬ 
sions  using  route  k  after  S(n+1)  missions  were  scheduled.  Then  the  problem  of 
picking  a  mission  type  to  minimize  expression  (A2),  i.e., 

min(Z(n+l,l) , . . .  ,Z(n+l  ,R) ) , 

is  equivalent  (as  far  as  picking  the  route)  to 

min(e(n+l ,1) ,...,  e(n+l,R)).  (A3) 


Remarks: 


1.  The  solution  in  general  is  not  unique. 

2.  (A3)  is  equivalent  to 

max( |e(n+l ,k) J :  e(n+l,k)<0)  (A4) 

since  e(n+l,l)+  ...  +  e(n+l,R) 

=  (S(n,l)-p(l)S(n+l) )  +  ...  +  (S(n ,R)-p(R)S(n+l) ) 

=  ( S(n , 1 )-p( 1 )S(n )-p( 1 ) )  +  ...  +  (S(n,R)-p(R)S(n)-p(R)) 

=  -(p(I  V-.-+P(R)) 

=  -1, 

so  e(n+l,k)<0  for  at  least  one  k. 

3.  In  terms  of  (A4),  the  conclusion  of  the  lemma  is  quite  intuitive.  As 
e(n+l,k)  is  a  measure  of  the  difference  between  the  actual  and  the  expected 
number  of  type-k  missions  scheduled,  the  lemma  says,  first,  that  one  should 
ignore  those  e(n+l,k)  that  are  positive;  this  makes  sense,  since  e(n+l,k)>0 
implies  that  type-k  missions  have  been  overscheduled.  Secondly,  it  says  that 
among  underscheduled  missions  (e(n  +  l,k)<0) ,  the  one  most  underscheduled  should 
be  chosen. 

Proof  of  Lemma  A.  If  a  type-k  mission  is  chosen,  then  the  average  absolute 
deviation  from  S(n+1)  is 

Z(n+l,k)  =  p(l) j(S(n,l)+v(l,k))/p(l)  -  S(n+l)|+  ... 

+p(R)|(S(n,R)+v(R,k))/p(R)  -  S(n+l)|. 
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Now  Z(n*l,k)  can  also  be  written  as 


Z(n+1 ,k)  =  p ( 1 ) |S(n,l)/p(l)  -  S(n+1)  +  ...  +  p(R)|S(n,R)/p(R) 

-  S(n+1)  +  p(k)  (S(n;k)+l)/p(k)-S(n+l) 

-  p(k) |s(n,k)/p(k)  -  S (n+1 ) j . 

Ignoring  the  terms  that  do  not  depend  on  k  (and  thus  do  not  contribute  any¬ 
thing  to  the  minimization  process),  the  original  minimization  is  equivalent  to 
(as  far  as  picking  the  route) 

min( |e(n+l,k)  +  lj  -  |e(n+l,k)|)  (A5) 

where 

e(n+l,k)  =  S(n,k)  -  p(k)S(n+l). 

Now  the  function  f(x)  =  jx  +  lj  -  jxj  is  nondecreasing  (-1  if  x  <  -1;  2x+l 
if  x  is  between  -1  and  0;  and  1  if  x  >  C);  hence  (A5)  is  equivalent  to 
mi n(e(n+l , 1 ) ,e(n+l ,2) , . . . ,e(n+l,R ) ) . 

Minimization  of  Average  Squared  Deviation 

Instead  of  minimizing  (A2),  here  we  minimize  p(I)  x  (S(n,l)/p(l)-S(n))2  + 
—  +  p(r)  x  (S(n,R)/p(r)-S(n) ) - .  By  a  similar  argument,  we  can  show  that  the 
type-k  mission  to  be  selected  at  time  n+1  (if  any)  is  one  for  which 

(2S(n, j )+l )pj  (j=l,...,R) 

is  minimum.  We  have  used  this  criterion  in  our  sample  problem. 
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APPENDIX  B:  NOMENCLATURE 


A 

B 

B(i ) 
B(l,n) 
B(i ,n,k , 

C 

c  ( i ) 

C(i  ,n) 
D(l,n) 
D(i ,n,k, 

d(i ,n) 
e(n,k) 

F (i ,n,k, 

f (n,i ,k) 

G 

l(i) 

M 

M(n) 

M(i ,n) 

m 

N 

i> 

P(i) 


The  family  of  arcs  in  the  graph,  G,  representing  the  airlift 
system. 

The  number  of  airbases  in  the  system. 

The  i-th  airbase  in  the  system.  (B(i)  is  the  home  base.) 

A  route  node  representing  the  home  base  at  time  n. 

h)  A  route  node  representing  airbase  B(i)  at  time  n;  this  is  the  h-th 
stop  on  the  k-th  route. 

The  number  of  aircrews  in  the  system. 

The  number  of  aircrews  initially  staged  at  airbase  B(i). 

An  aircrew  node  representing  airbase  B(i)  at  time  n. 

The  number  of  aircraft  delayed  at  home  base  at  time  n. 

h)  The  number  of  aircraft  delayed  at  airbase  B(i)  at  time  n  and  on 

the  h-th  stop  of  the  k-th  route. 

The  number  of  aircrews  delayed  at  airbase  B(i)  at  time  n. 

The  difference  between  the  actual  and  expected  number  of 

missions  of  type  k  scheduled  at  time  n. 

h)  The  number  of  aircraft  leaving  airbase  B(i)  at  time  n,  where 
B(i)  is  the  h-th  stop  on  the  k-th  route. 

The  number  of  aircrews  leaving  B(i)  towards  airbase  B(k)  at 
time  n. 

The  graph  representing  the  airlift  system. 

The  number  of  legs  in  the  i-th  route.  There  are  L(i)  +  1  stops 
in  the  i-th  route,  B(l)  being  the  first  and  the  last. 

The  master  schedule  node. 

The  number  of  aircraft  starting  maintenance  at  time  n. 

A  subsidiary  schedule  node  related  to  the  i-th  route  at  time  n. 

Aircraft  maintenance  time. 

The  family  of  nodes  in  the  graph,  G,  representing  the  system. 

The  number  of  aircraft  in  the  system. 

The  proportional  usage  of  route  i. 


I? 


R  The  number  of  routes  in  the  system. 

R(n)  A  set  of  route  nodes  related  to  airbases  at  time  n. 

r  The  aircrew  rest  time. 

r(i,n)  The  number  of  aircrew  resting  at  airbase  B(i)  at  time  n. 

S  The  source  node  in  the  graph,  G,  representing  the  system. 

S(n)  The  total  number  of  missions  sent  out  as  of  time  n. 

S(n,k)  The  total  number  of  type-k  missions  that  were  sent  ouc  as  of 

time  n. 

s(n)  The  total  number  of  missions  that  are  sent  out  at  time  n. 

s(n,k)  The  total  number  of  type-k  missions  that  are  sent  out  at 

time  n. 

T  Duration  of  airlift  being  considered. 

t(i,j)  Flight  time  between  airbases  B(i)  and  B(j). 

1)  Aircraft  utilization;  i.e.,  the  total  number  of  flying  hours  by 

all  aircraft  during  time  period  T. 

v(i,j)  A  number  equal  to  1  when  i  =  j  and  to  0  otherwise. 

w(k,h)  Flight  time  between  airbases  b(k,h)  and  b(k,h+l). 
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